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- Abstract: We study the sources of systematic errors in the measurement of the Z 

psj . cross-sections at the LHC. We consider the systematic errors in both the total cross-section 

i and acceptance for anticipated experimental cuts. We include the best available analysis 

g ; of QCD effects at NNLO in assessing the effect of higher order corrections and PDF and 

QQ I scale uncertainties on the theoretical acceptance. In addition, we evaluate the error due to 

' missing NLO electroweak corrections and propose which MC generators and computational 
schemes should be implemented to best simulate the events. 
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1. Introduction 

A precise measurement of gauge boson production cross-sections for pp scattering will be 
crucial at the LHC. W and Z bosons will be produced copiously, and a careful measurement 
of their production cross-sections will be important in testing the Standard Model more 
rigorously than ever before, and uncovering signs of new physics which may appear through 
radiative corrections. In addition, these cross-sections have been proposed as a "standard 
candle" for measuring the luminosity through a comparison of the measured rates to the 
best theoretical calculations of the cross-section. Investigation of this means of measuring 
luminosity began at the Tevatron and will continue at the LHC |Q, ^. 

In this paper, we will concentrate on Z production. A comparable analysis of sys- 
tematic uncertainties in W production appeared in Ref. @. Since that time, both the 
experimental approach to the measurements and the theoretical results needed to calcu- 
late them have both been refined. NNLO QCD calculations of these processes, previously 
available only for the total cross-section and rapidity distribution H], are now available 
in differential form permitting an analysis of the effect of experimental cuts on the 
pseudorapidity and transverse momentum of the final state leptons. 

The high luminosity (10^^ - 10^'' cm^s"^) at the LHC insures that systematic errors 
will play a dominant role in determining the accuracy of the cross-section. Thus, we present 
an analysis of the effect of the theoretical uncertainty in the evaluation of the acceptance, 
and propose which among the various available MC generators and computational schemes 
should be implemented to best simulate the events. 

This paper is organized as follows. Sec. § will give an overview of the calculation and 
the computational tools used in the analysis. The next four sections are each devoted to 
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estimating a class of systematic errors: electroweak corrections in Sec. ^, NNLO QCD in 
Sec. QCD scale dependence in Sec. ^, and parton distribution function uncertainties in 
Sec. ^. Finally, the results are compiled and summarized in Sec. ^ 

2. Theoretical Calculations and MC Generators 

The dominant production mechanism for Z bosons is the Drell-Yan process 0, in which 
a quark and antiquark annihilate to form a vector boson, which subsequently decays into 
a lepton pair. It is through this lepton pair that the production process is observed. 
Such pairs may as well be produced via an intermediate photon 7*, so both cases should be 
considered together. The Z production cross-section may be inferred experimentally from 
the number N'i^^ ^ of observed events via the relation 

AT^bs^^ ^ ^tot ^Yi{Zh* ^ tl-)A2i^* j Cdt. (2.1) 

acceptance obtained after applying the experimental selection criteria. For 
example, if the cuts require > p™*", < r/^ < rj^^^, then 

1 rV~s/'i rVmax + 

dp^* dpr^'^ dr](^+dr]£- 

In practice, further cuts on the invariant mass M^g of the lepton pair may be included 
to prevent the cross-section from being dominated by photons, which give a divergent 
contribution at low energies. 

Alternatively to the Z production cross-section measurement, the corrected Z yield 
can be used as a standard candle for a luminosity monitor in LHC if one calculates the 
cross-section and solves for J Cdt. The theoretical cross-section may be constructed by 
convoluting a parton-level cross-section aab for partons a and b with the parton density 
functions (PDFs) fa, ft for these partons, 

a'^'ipp ^ Z/7* ^ £+£-) = Y, [' dxidx2fa{xi)Mx2) aab{xi,X2), (2.3) 



a,b 







integrating over the momentum fractions xi,X2, and applying cuts relevant to the exper- 
iment. Theoretical errors come from limitations in the order of the calculation of aab, 
on its completeness (for example, on whether it includes electroweak corrections or 'j* /Z 
interference, and on whether any phase space variables or spins have been averaged), and 
from errors in the PDFs. ^ 



^The representation eq. of the cross-section is intended to illustrate the manner in which the mea- 
surement may be used to infer the Z production cross-section, and does not imply that a narrow resonance 



approximation is actually used in calculating the cross-section 2.c 
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Since the final state may include additional partons which form a shower, the output 
from the hard QCD process must be fed to a shower generator to generate a realistic final 
state seen in a detector. This is possible only if the cross-section is simulated in an event 
generator. Calculating the acceptance for all but the simplest cuts will normally require 
an event generator as well. 

Thus, when constructing a simulation of an experiment, there is a range of choices 
which can be made among the tools currently available. An efficient calculation requires 
selecting those adequate to meet the anticipated precision requirements, without perform- 
ing unnecessarily complex calculations. For example, while NNLO calculations are now 
available, the cross-sections are very complicated, do not always converge well, and require 
substantial time to calculate. For certain choices of cuts, it may be found that the effect of 
the NNLO result can be minimized, or that it can be represented by a simplified function 
for the parameters of interest. We will compare several possible schemes for calculating 
the Z production cross-section and acceptance, and consider the systematic errors arising 
for these schemes. 

The most basic way to generate events is through one of the showering programs, such 
as PYTHIA @, HERWIG [|, IS A JET [|lO| or SHERPA These vary somewhat in 

their assumptions and range of effects included, but they all start with hard partons at a 
high energy scale and branch to form partons at lower scales, which permits a description 
of hadronization and realistic events. On their own, these programs typically rely on a 
leading order hard matrix element and include only a leading-log resummation of soft 
and collinear radiation in the shower, limiting their value in describing events with large 
transverse momentum. In addition, ISAJET lacks color-coherence, which is important in 
predicting the correct distribution of soft jets fl^ . 

Fully exclusive NLO QCD calculations are available for W and Z boson production |l3| . 
The MC generator MC@NLO @ combines a parton-level NLO QCD calculation with the 
HERWIG parton shower, thus removing some of the limitations of a showering program 
alone. 

Since a ~ at LHC energies, NLO electroweak (EWK) corrections should appear 
at the same order as NNLO QCD. The MC@NLO package is missing EWK corrections, 
but the most important of these corrections under the Z peak is expected to be QED final 
state radiation 0. This can be obtained by combining MC@NLO with PHOTOS |jl5|, an 
add-on program which generates multi-photon emission from events created by the host 
program. Another program, HORACE [|l8|, is available which includes exact 0{a) EWK 
corrections together with a final state QED parton shower. 

The other available NLO and NNLO calculations are implemented as MC integrations, 
which can calculate a cross-section but do not provide unweighted events. Some of these 
are more differential than others. For example, the NNLO rapidity distribution is available 
in a program Vrap |^, but this distribution alone is not sufficient to calculate acceptances 
with cuts on the lepton pseudorapidities and transverse momenta. A differential version 
of this NNLO calculation is implemented in a program FEWZ but this is not an event 
generator. Another available program is ResBos-A |1£], which resums soft and collinear 
initial state QCD radiation to all orders and includes NLO final state QED radiation. 
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Resummation gives this program an advantage in realistically describing the small 
regime. 

Our analysis is conducted for di-lepton final states. The available calculations typically 
set the lepton masses to zero, so the lepton masses will be neglected throughout this paper 
and the choice of final state lepton has no effect on the calculations. In all results, i may 
be interpreted as either an electron or muon. We have chosen three sets of experimental 
cuts to reflect detector capabilities and to demonstrate the impact of physics effects on the 
acceptances depending on the selection criteria. 



3. Electroweak Corrections 

As noted above, both NNLO QCD corrections and NLO electroweak (EWK) corrections 
are expected to be needed to reach precisions on the order of 1% or better in Z boson 
production. NLO electroweak [pO|] and QCD corrections |jl3| are known both for W boson 
production and Z boson production. However, current state-of-the-art MC generators do 
not include both sets of corrections. The generator MC@NLO ||l^ combines a MC event 
generator with NLO calculations of rates for QCD processes and uses the HERWIG event 
generator for the parton showering, but it does not include EWK corrections. Final state 
QED can be added using PHOTOS a process-independent module for adding 

multi-photon emission to events created by a host generator. However, some 0(a) EWK 
corrections are still missing. 



To study the error arising from missing 0{a) EWK corrections, we used HORACE |18|, 
a MC event generator that includes initial and final-state QED radiation in a photon shower 
approximation and exact 0{a) EWK corrections matched to a leading-log QED shower. 
To determine the magnitude of the error, we then compared the results from this generator 
to a Born- level calculation with final-state QED corrections added by PHOTOS. 

Specifically, we compared pp — > Z/j* — > events generated by HORACE with 

the full 0(a) corrections and parton-showered with HERWIG, to these events generated 
again by HORACE, but without EWK corrections (Born-level), showered with HER- 
WIG+PHOTOS. CTEQ6.5M parton distribution functions [^T| were used in the calcu- 
lations. The results are shown in Table |^ and in Figs. |l| ^ 

For Table the events are generated in the kinematic region defined by final state 
invariant mass M« > 40 GeV/c^, and for each lepton, p^ > 5 GeV/c, and \7]^\ < 50.0. 
Two sets of cuts are used: 

Loose Cut: M« > 40 GeV/c^ > 5 GeV/c, |r?J < 50.0, 

Tight Cut: 40 < M« < 140 GeV/c^, > 20 GeV/c, |r?J < 2.0. 

The first row in the table shows the total generator-level cross-sections before QCD 
parton showering (identified by the label "No PS"). At Born level, without the corrections, 
the "loose" cut is essentially identical to the generator-level cut. The Born-I-FSR column 
shows the effect of applying final state radiation (FSR) corrections only via PHOTOS. 



In PHOTOS, FSR affects the rates through the cuts only. As noted in Ref. |16], the 



combined effect of the complete first order real and virtual corrections not included in the 
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Photonic and Electroweak Corrections 





Born 


Born+FSR 


ElectroWeak 


Difference 


a (No PS) 
a (Loose Cut) 
cj(Tight Cut) 
A (Loose Cut) 
A (Tight Cut) 


1984.2 ± 2.0 
1984.2 ± 2.0 
612.5 ± 1.1 
0.9999 ± 0.0000 
0.3087 ± 0.0005 


1984.2 ± 2.0 
1964.6 ± 2.0 
597.6 ± 1.1 
0.9901 ±0.0001 
0.3012 ±0.0005 


1995.7 ±2.0 
1961.4 ±2.0 
595.3 ± 1.1 
0.9828 ± 0.0001 
0.2983 ± 0.0005 


0.58 ± 0.14% 
0.16 ± 0.14% 
0.38 ± 0.26% 
0.74 ± 0.02% 
0.96 ± 0.21% 



Table 1: Calculation of the Z/^* I'^l cross-section a and acceptance A for various EWK 
corrections generated using HORACE 3.1, for £ = e or ^. 



Generated Z -> Z+Z^lny) Mass Generated l*r Mass 




40 50 60 70 80 90 100 110 120 130 140 40 50 60 70 80 90 100 110 120 130 140 



{GeV/c=^) M„ (GeV/c=^) 

(a) (b) 

Figure 1: Comparison of the (a) Z boson invariant mass distributions and (b) i'^i" invariant mass 
distributions for the process Z/j* — > £^i~{n^) in HORACE 3.1 including 0{a) EWK corrections 
showered with HERWIG (open red squares), HORACE Born-level showered with HERWIG plus 
PHOTOS (black circles), and HORACE Born-level (blue stars). 



standard version of PHOTOS would increase the total rate by a factor of 1 ± 3a/47r, a 
0.17% increase. The ElectroWeak column includes the full HORACE EWK corrections. In 
the final column, we give the difference between the previous two columns, to compare the 
full EWK correction FSR alone. The results show agreement within 1% between the two 
schemes. The maximum error in the cross-section is 0.58% without added parton shower, 
and the maximum error in the acceptance is 0.96% for the tighter cut. Therefore, we 
recommend using MC@NLO interfaced with PHOTOS as our primary event generator for 
measurements at the Z peak, until higher precision is required [p^ . 

For completeness, we also compared MC@NLO interfaced with PHOTOS with dis- 
tributions from ResBos-A |]l9[, a MC simulation that includes final state NLO QED cor- 
rections to W/Z boson production and higher-order logarithmic resummation of soft and 
collinear QCD radiation. Fig. ^ shows the distributions for and the of the Z boson, 
respectively. The latter exhibits especially well the effects of the soft and collinear resum- 
mation at low in ResBos-A [19|. The only generator- level cut used in these comparison 
plots is the Mu > 50 GeV/c^. Otherwise, the default parameters were used with no tuning. 
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Number of FSR Photons 



FSR Photon Transverse Momentum 



HOR Born+Herwig 
HOR Born+Herwig+PHOTOS 
HOR EWK+Herwig 



> 




10 



Figure 2: Comparison of the number n 
of final state radiation (FSR) photons in 
Z/7* ^ e+i-{n-/) for HORACE 3.1 in- 
cluding 0{a) EWK corrections showered 
with HERWIG (open red squares), HO- 
RACE Born-level showered with HERWIG 
plus PHOTOS (black circles), and HORACE 
Born-level (blue stars). 



HOR Born+Herwig+PHOTOS 
HOR EWK+Herwig 




2 4 6 8 10 12 14 

Pt, (GeV/c) 

Figure 3: Comparison of Z/j* ^+^"(717) 
final state radiation (FSR) transverse mo- 
mentum distributions for HORACE 3.1 in- 
cluding EWK corrections showered with 
HERWIG (open red squares) and HORACE 
Born-level showered with HERWIG plus 
PHOTOS (black circles). 



Invariant t r Mass Generated Z Transverse Momentum 




M„-(GeV/c2) p^(GeV/c) 
(a) (b) 

Figure 4: Comparison of the (a) di-lepton invariant mass and (b) Z transverse momentum distri- 
butions for the process Z/j* i+i-{n-/), in MC@NLO with PHOTOS (red squares) and Resbos-A 
(black circles). 



4. NNLO QCD Uncertainties 

QCD uncertainties include errors due to missing higher-order corrections in the hard matrix 
element, uncertainties in the parton distribution functions, and approximations made in the 
showering algorithms. In the following, we will evaluate the errors introduced by omitting 
the NNLO corrections by using MC@NLO, and calculate iC-factors which can be used 
to introduce NNLO corrections to the MC@NLO calculation. We will also examine the 
effect of uncertainties in the PDFs. For these studies, we choose three sets of experimental 
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Invariant Mass (GeV/c^) 


Pseudorapidity 


Transverse Momentum (GeV/c) 


Cut 1 


Mu > 40 


|r?J < 2.0 


> 20 


Cut 2 


Mu > 40 


1.5 < IryJ < 2.3 


>20 


Cut 3 


79 < M« < 104 


Kl <2.0 


>20 



Table 2: Acceptance regions for final state leptons in NNLO studies 



cuts 1^ to reflect detector capabilities and to demonstrate the impact of physics effects on the 
acceptances depending on the selection criteria. Here, r/^ andp^ are the pseudorapidity and 
transverse momentum of the final state leptons. The different rapidity ranges and invariant 
di-lepton mass cuts provide useful separation for between regions of the Z spectrum which 
have different sensitivities to some of the sources of uncertainties and evaluate the impact 
of mass cut on the theoretical error. 

We begin by examining the NNLO corrections, using the state-of-the-art program 
FEWZ IQ, which is differential in the di-lepton invariant mass, and the lepton transverse 
momenta and pseudorapidities. The FEWZ program is at NNLO in perturbative QCD, and 
fully differential, giving correct acceptances including spin correlations, as well as taking 
into account finite widths effects and j* — Z interference. Since we are interested primarily 
in studies about the Z peak, we choose the renormalization and factorization scales to be 
fj,^ = fj,^ = Mz- Scale dependence will be discussed in detail the next section. 

A comparison of the effects of higher order QCD corrections on the cross-section and 
acceptance is presented in Table ^ and Figs ^ - Both the NLO and NNLO calculations 
are done with CTEQ6.5M PDFs |21|. We comment further on this choice at the end of 
this section. All results in the table are calculated at scale Mz- In the figures, the NLO 
results are displayed as a band spanning the range of scales from AIz/2 to 2Mz- The scale 
dependence of the NNLO result is small enough to be comparable to the precision of the 



NNLO Cross Sections a (pb) 



Cut 


MC@NLO 


FEWZ NLO 


FEWZ NNLO 


ii'-factor 


Kx MC@NLO 


^tot 


2331 ± 3 


2358.1 ±2.3 


2334.9 ±4.6 


0.9902 ± 0.0022 


2308 ± 6 


1 


703.6 ± 1.1 


716.0 ±0.7 


726.2 ± 5.5 


1.0142 ±0.0077 


713.6 ±5.6 


2 


71.3 ±0.3 


74.1 ± 0.07 


73.39 ± 1.96 


0.9902 ± 0.0280 


70.6 ± 1.9 


3 


623.5 ± 1.0 


657.2 ±0.7 


650.4 ± 4.0 


0.9897 ±0.0062 


617.1 ±4.0 






NNLO Acceptances 


;%) 




Cut 


MC@NLO 


FEWZ NLO 


FEWZ NNLO 


A'- fact or 


Kx MC@NLO 


1 


30.18 ±0.06 


30.36 ± 0.04 


31.10 ±0.24 


1.0243 ± 0.0081 


30.92 ± 0.25 


2 


3.06 ±0.01 


3.143 ±0.004 


3.143 ± 0.084 


1.0000 ± 0.0268 


3.06 ± 0.08 


3 


26.75 ± 0.06 


27.87 ±0.04 


27.86 ±0.18 


0.9995 ± 0.0066 


26.73 ±0.19 



Table 3: Calculation of the Z/j* £+1^ {£ = e or fi) cross-section at NLO using MC@NLO, 
and at NLO and NNLO using FEWZ, for the cut region defined in Table The acceptances are 
relative to a total generated cross-section a'°* for Mu > 40 GeV/c^ . 
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Lepton Cut: (M^^> 40, |ri| < 2.0) 




Lepton p^ Cut: (M^^> 40, (nl < 2.0) 

0.35 hi- ' ; T 




15 20 25 30 

, (GeV/c) 
(b) 



Figure 5: Cross-section (a) and acceptance (b) versus cut on lepton at NLO (hashed bands) 
and NNLO (points), as calculated using FEWZ. 



Lepton |ri| Cut: (M^^> 40, p.^ > 20) 



Lepton |ti| Cut: (M^^> 40, p^ > 20) 




(a) 



ra 




(b) 



Figure 6: Cross-section (a) and acceptance (b) versus cut on lepton \r]\ at NLO (hashed bands) 
and NNLO (points), as calculated using FEWZ. 



Invariant Mass Cut: (p.^ > 20, tnl < 2.0) 



Invariant Mass Cut: (p^ > 20, (n| < 2.0) 




AM„(GeV/c^) 

(a) 




AM„(GeV/c^) 
(b) 



Figure 7: Cross-section (a) and acceptance (b) versus cut on di-lepton invariant mass Mu, at 
NLO (hashed bands) and NNLO (points), as calculated using FEWZ. M« is constrained to be in 
a band of width AMu centered at Mz- 
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MC evaluation of the integrals, so only the average of the high and low scales is plotted, 
with error bars reflecting a combination of statistical and scale variation uncertainties. 
Fig. ^ shows the cross-section and acceptance for an invariant mass band of width AM^^ 
centered at Mz-, to study the effect of selecting different cuts about the Z peak. 

Since the NNLO matrix element has not yet been interfaced to a shower, we cannot 
directly compare FEWZ to MC@NLO. The best we can do at this time is to use MC@NLO 
to obtain the NLO showered result, and multiply this by a X-factor obtained by taking 
the ratio of the NNLO to NLO results derived from FEWZ. This procedure is reasonable 
except in threshold regimes where the fixed-order NLO result in FEWZ is unreliable. In 
fact, the NLO cross-section calculated by FEWZ can become negative near the threshold 

= Mz/2. Similar methods have been used for calculating NNLO corrections to Higgs 
production |2^] . The differences of these /C- factors from unity are shown in Fig. ^ for both 
the cross-sections and acceptances, as a function of cuts on the lepton and rj as well 
as on the width AM^^ of an invariant mass cut centered at Mz- The resulting accepted 
cross-section is shown in the Ky. MC@NLO column of Table ^ and in Fig. ^ as a function 
of the same cuts as in Fig. |. The size of -fC — 1 is a good indicator of the error due to 
missing NNLO if MC@NLO is used without corrections. 

The results in table ^ show X-factors corresponding to an NNLO correction of about 
1% for the cross-sections, or up to 2.4% for the Cut 1 acceptance. FEWZ's convergence 
was not as good for Cut 2, however, leading to an almost 3% technical error in K from 
the evaluation.^ In the other cases, a 1% evaluation of the NNLO correction was generally 
possible within a reasonable run time, which was typically of order one month on the 
clusters used in this study. Longer run times do not appear to improve the convergence 
significantly. The convergence limitations are due to the high dimensionality (11) of the 
Vegas |2^ integrals and strong peaking of the integrands. 

We should note that CTEQ6.5M structure functions were used in the above calculations 
because the K factors calculated are intended for rescaling an MC@NLO calculation using 
CTEQ PDFs. However, since these PDFs do not include NNLO corrections, it is useful 
to note the effect of repeating the calculations with MRST 2002 NNLO PDFs. The K- 
factors for the cross-section and acceptance for each of the three cuts is shown in Table 
^ for both sets of PDFs, for comparison. The i^(c'")-factors are ratios of NNLO to NLO 
cross-sections, and K{A) are ratios of NNLO to NLO acceptances relative to a total cross 
section for > 40 GeV. Except for the cross-section K for Cut 1, where agreement is to 
1.8%, the two PDF choices agree to within 1%, and within the computational uncertainty. 



^Here, "technical error" refers to a limitation on the precision arising from the specific calculation of 
the process, rather than from the approximations made in the calculation at a theoretical level (NNLO vs. 
NLO, for example). Technical precision is discussed further in the conclusions. 
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Cross-Sections 



Acceptances 



0.05 

0.04 
0.03 
0.02 
0.01 


-0.01 
-0.02 
-0.03 
-0.04 
-0.05, 



Lepton Cut: (M^^> 40, |ti| < 2.0) 



Pt, 



20 25 

:„ (GeV/c) 



O.lp- 
0.08- 
0.06 
0.04 H 
0.02 



< 



(a) 



Lepton p^ Cut: (M^^> 40, Ml < 2.0) 



711 1 1 1 1 mi ^ I 



-0.02 r 
-0.04 ~- 
■0.06 
■0.08 



Pt, 



20 25 

, (GeV/c) 



Lepton |ti| Cut: (M^^> 40, p^ > 20) 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 



0.1 

0.08 - 

0.06 

0.04 H 

0.02 ^ 
0^ 
-0.02 r- 
-0.04r- 
-0.06 r- 
-0.08 r- 

-O.f 



Lepton |ri| Cut: (M^^> 40, p^ > 20) 



2.5 

Inl 

I II „ 



(b) 
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-0.05 



Invariant Mass Cut: (p^ > 20, hi < 2.0) 



I , ■ 






I 














1- 




















ii 



































92 94 96 98 100 102 



AM„{GeV/c'') 



0.05 p- 
0.04 H - 
0.03 ^ 
0.02^ 
0.01 E 



< 



(c) 



Invariant Mass Cut: (p^ > 20, |ri| < 2.0) 



-0.01^ 
■0.02 H 
■0.03 H 
■0.04^ 
■0.05^ 



96 98 

AM„(GeV/c^) 



Figure 8: Fractional difference in the NNLO and NLO cross-sections (left-hand side) and accep- 
tances (right-hand side) as a fmiction of the lepton (a) p^, (b) \ri\, and (c) M« cuts as in Figs. 
Figs 1^-0- These differences are the factor K — 1 for the cross-section and acceptances, respectively. 
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Lepton p Cut: (M > 40, |ri| < 2.0) 



800 

700 
600 
500 
400 
300 
200 



j: X z i 



Lepton |ri| Cut: (M^^> 40, p^ > 20) 
1600 r 

1200^ 



J3 

Q. : 
"-i 1000- 



600: 
400: 
200 : 



10 15 20 25 30 35 40 

p^^JGeV/c) 
(a) 



0.5 1 1.5 2 2.5 3 3.5 

hi 

' " max 

(b) 





4 4.5 5 



Invariant Mass Cut: (p > 20, |ri| < 2.0) 



700 - 
650 - 



S~ 600 - 
Q. 



500 r 
450 r 
400 EL 



, I I ■ 



94 96 98 100 

AM„(GeV/c^) 



(c) 



Figure 9: Accepted NNLO cross-section estimated from MC@NLO scaled by the if-factor versus 
the lepton (a) p^, (b) |?7|, and (c) M« cuts as in Figs. ||-^- 



K-factors with CTEQ and MRST PDFs 



Cut 


CTEQ K{a) 


MRST K{a) 


CTEQ K{A) 


MRST K{A) 


1 

2 
3 


1.0142 ±0.0077 
0.9902 ± 0.0280 
0.9897 ± 0.0062 


0.9963 ± 0.0109 
0.9830 ± 0.0348 
0.9851 ± 0.0082 


1.0243 ±0.0081 
1.0000 ±0.0268 
0.9995 ± 0.0066 


1.0211 ±0.0113 
1.0074 ±0.0357 
1.0096 ±0.0087 



Table 4: Comparison of the K factors for cross-sections and acceptances calculated using 
CTEQ6.5M (NLO) and MRST 2002 (NNLO) PDFs. K{a) is the isT-factor for the cross-section 
and K{A) is the iiT-factor for the acceptance, as displayed in Table ^. 



5. Scale Dependence 

Perturbative QCD calculations at fixed order depend on the factorization and renormal- 
ization scales introduced in the calculation. Thus, the previous calculations have an added 
uncertainty due to the choice of certain fictitious scales appearing in the calculation. In a 
complete, all order calculation, there would be no dependence on these scales. However, 
in a fixed order calculation matched to PDFs, a dependence on the factorization scale 
and renormalization scale appear in the final results. The effect of scale choice is 
significant at NLO. Adding NNLO effects is found to reduce scale dependence consider- 
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ably H, g, though it can remain significant near thresholds where large logarithms render 
a fixed-order result unreliable. 

As is customary, we will choose the renormalization and factorization scales to be 
identical, and investigate the scale dependence by varying them by a factor of 2 or 1/2 
about a central value of fJ-p = Mz, which is typical of the scales in our acceptance, and 
was the central value chosen in the previous section. 

Table ^ included only the central scale Mz- Tables ^ and ^ and show the total cross 
sections and acceptances for lepton production calculated by FEWZ at three different 
renormalization and factorization scales Mz/2, Mz, and 2Mz- The acceptances for the 
final state leptons are as defined in Table ^. For a measure of the size of the scale de- 
pendence, the final column of each table shows the maximum difference between the three 
values divided by average, with an error calculated assuming the statistical errors in the 
three MC runs in each row are independent. 

We can see that the scale dependence of the cross-sections at NLO is typically of 
order ±2.5%. The scale dependence of NLO acceptances is dramatically reduced due to 
correlations in the scale dependence of the cut and uncut cross-sections used to compute 
it. Adding NNLO reduces the scale dependence of the cross-sections to ±1.2% or less. 
These estimates were obtained using the CTEQ6.5M NLO PDFs, but are compatible with 
the scale dependence found using MRST 2002 NNLO PDFs for the same calculations, 
within the limits of the computational errors. A better evaluation of the NNLO result 
may show the scale dependence to be even smaller, since it is at the same level as the MC 
precision. This is particularly true for Cut 2, for large values of \rj\, where the cross-section 
is relatively small and the MC errors relatively large. Improving the convergence of FEWZ 
would reduce these errors. 



6. Uncertainties Due to the Parton Distribution Function 



Phenomenological parameterizations of the PDFs are taken from a global fit to data. There- 
fore, uncertainties on the PDFs arising from diverse experimental and theoretical sources 
will propagate from the global analysis into the predictions for the W/Z cross-sections. 
Figure 10 shows the results of the inclusive Z to di-lepton production cross-section us- 
ing various CTEQ H and MRST Q PDFs. The upward shift of about 7% (between 
CTEQ6.1 and 6.5 and MRST2004 and 2006) results from the inclusion of heavy quark 
effects in the latest PDF calculations. The acceptance due to the cuts in Table ^ using 



each of these PDFs is shown in Fig. 11 



The uncertainties in the PDFs arising from the experimental statistical and systematic 
uncertainties, and the effect on the production cross-section of the Z boson, have been 
studied using the standard methods proposed in Refs. [21, 24]. For the standard set of 



^Some theoretical issues which may affect the contribution of the PDFs to the NNLO A"-factor are 
not included here as we are concerned primarily with the error at NLO. See Refs. |^] for details. The 
comparison at the end of Sect. ^ suggests that such effects are likely to be small in the K factors calculated 
here. 
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Total Cross-Section (in pb), M« > 40 GeV 



Order 


Mz/2 Mz 2Mz 


Aa/a 


NLO 
NNLO 


2297.4 ± 2.3 2358.1 ±2.3 2418.3 ± 2.3 
2320.4 ±4.6 2334.9 ±4.6 2350.1 ±4.6 


0.0513 ±0.0014 
0.0127 ± 0.0028 


Cut Region 1 


Order 


Mz/2 Mz 2Mz 


Act/ct 


NLO 
NNLO 


698.2 ±0.7 716.0 ±0.7 735.2 ± 0.7 
718.9 ±6.1 726.2 ±5.5 719.8 ±4.8 


0.0516 ±0.0014 
0.0101 ±0.0108 


Cut Region 2 


Order 


Mz/2 Mz 2Mz 


Aa/W 


NLO 
NNLO 


72.55 ± 0.07 74.12 ± 0.07 75.78 ± 0.08 
73.78 ± 2.00 73.39 ± 1.96 75.10 ± 1.52 


0.0436 ± 0.0014 
0.0231 ± 0.0351 


Cut Region 3 


Order 


Mz/2 Mz 2Mz 


Aa/W 


NLO 

NNLO 


640.2 ±0.6 657.2 ±0.7 673.9 ±0.7 

661.3 ±4.7 650.4 ±4.0 656.2 ± 3.8 


0.0513 ±0.0014 
0.0166 ± 0.0090 



Table 5: Scale dependence of the total and accepted cross-sections (in pb) for Z boson production 
with Mu > 40 GeV/c^ calculated by the FEWZ program at order NLO and NNLO. The final 
column is a measure of scale dependence obtained by dividing the maximum spread by the average 
for the three points. 





Cut Region 1 (% 


Accepted) 




Order 


Mz/2 Mz 


2Mz 


AA/A 


NLO 


30.39 ± 0.04 30.36 ± 0.04 


30.40 ± 0.04 


0.0013 ± 0.0020 


NNLO 


30.98 ± 0.27 31.10 ± 0.21 


30.40 ± 0.04 


0.0153 ±0.0111 




Cut Region 2 (% 


Accepted) 




Order 


Mz/2 Mz 


2Mz 


AA/A 


NLO 


3.158 ±0.031 3.143 ±0.030 


3.134 ±0.033 


0.0077 ± 0.0020 


NNLO 


3.180 ±0.086 3.143 ±0.084 


3.196 ±0.065 


0.0165 ±0.0353 




Cut Region 3 (% 


Accepted) 




Order 


Mz/2 Mz 


2Mz 


AA/A 


NLO 


27.87 ±0.04 27.87 ±0.04 


27.87 ± 0.04 


0.0001 ± 0.0020 


NNLO 


28.50 ±0.21 27.86 ±0.18 


27.92 ±0.17 


0.0229 ± 0.0095 



Table 6: Scale dependence of the acceptances A in the various cut regions for Z boson production 
with Mte > 40 GeV/c^ calculated by the FEWZ program at order NLO and NNLO. The final 
column is a measure of scale dependence obtained by dividing the maximum spread of the three 
preceding columns by their average. 
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Z -> tr NLO Cross-Section 
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Figure 10: Comparison of Z/7* cross-sections for Mu > 40 GeV/c for several recent 

PDF calculations. 



PDFs, corresponding to the minimum in the PDF parameter space, a complete set of 
eigenvector PDF sets, which characterize the region nearby the minimum and quantify its 
error, have been simultaneously calculated. From the minimum set and these "error" sets 
we calculate the best estimate and the uncertainty for the Z cross-section. We do this using 
the asymmetric Hessian error method, [26| where the cross-section results from the various 



eigenvector PDF sets have been combined according to the prescriptions found in |21, 26|. 
Fig. |l^ list the results for the different PDFs and Table ^ summarizes the results of the 
latest CTEQ and MRST PDF sets. These calculations were done using MC@NLO. The 
difference in the uncertainties (approximately a factor of two) between the results obtained 
from the CTEQ and MRST PDF error sets is due to different assumptions made by the 
groups while creating the eigenvector PDF sets. Note that the first column of Table ^ 
should agree with the CTEQ6.5 entries in Table 0, and do within the uncertainties quoted 
in Table ^ although different event samples were used. 

Finally we study the sensitivity of the kinematic acceptance calculations to the uncer- 



tainties affecting the PDF sets. Figs. |12| and |13| show the systematic error on the production 
cross-sections as a function of the |ry| cut and minimum lepton for variations on the 
three types of cuts in Table |2[ The fractional uncertainties, shown in in the same figures, 
demonstrate that the relative uncertainty in the cross-section is very flat as a function 
of the kinematic cuts, until the region of extreme cuts and low statistics in the MC are 
reached. The corresponding uncertainty on the acceptance as a function of the kinematic 



cuts is shown in Figs 14 and 15. These show a similar dependence to the cross-section 



uncertainties, though the fractional errors are smaller. 
7. Conclusions 

To evaluate the overall contribution from theoretical uncertainties to both the cross-section 
and acceptance calculations for the decay mode Z/^* £^£~ (£ = e or fi) at the LHC 
(and for Mgi > 40 GeV/c^) we add the uncertainties from each of the sources considered in 
the preceding sections. We compile the errors assuming that the calculation is done with 
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Mff 


> 40 GeV/c2 


Cut Region 1 


PDF Set 


a fob) 


Acr+ 


Aa_ 


a (pb) 


Acr+ 


Aa_ 


A 


AA+ 


AA_ 


CTEQ6.5 


2330 


103 


104 


703.6 


21.6 


26.7 


0.302 


0.004 


0.004 


MRST2006 


2333 


42 


40 


712.6 


13.6 


16.0 


0.305 


0.004 


0.004 


CTE06 1 


2155 


123 


109 


652.1 


30.2 


29.5 


0.303 


0.007 


0.005 


MRST2004 (NNLO) 


2193 


41 


45 


672.4 


12.6 


17.9 


0.302 


0.004 


0.004 


MRST2004 (NLO) 


2223 


42 


46 


662.8 


12.5 


17.9 


0.302 


0.004 


0.004 






> 40 GeV/c2 


Cut Region 2 


PDF Set 


(J (vh) 


A(T+ 


Act- 


(T (pb) 


A(T+ 


A(7_ 


A 


Ayl+ 


AA- 


CTEQ6.5 


2330 


103 


104 


71.1 


2.2 


3.4 


0.0305 


0.0001 


0.0005 


MRST2006 


2333 


42 


40 


72.0 


0.1 


2.4 


0.0308 


0.0000 


0.0010 


CTE06 1 


2155 


123 


109 


66.3 


3.0 


4.1 


0.0307 


0.0006 


0.0014 


MRST2004 rNNLO"! 


2193 


41 


45 


69.0 


0.3 


2.1 


0.0310 


0.0002 


0.0007 


MRST2004 (NLO) 


2223 


42 


46 


68.7 


0.3 


2.0 


0.0313 


0.0002 


0.0007 




Mu 


> 40 GeV/c2 


Cut Region 3 


PDF Set 


a (pb) 


A(T+ 


Act- 


a (pb) 


Acr+ 


A(7_ 


A 


AA+ 


AA_ 


CTEQ6.5 


2330 


103 


104 


623.5 


20.8 


22.0 


0.268 


0.005 


0.003 


MRST2006 


2333 


42 


40 


634.1 


10.4 


15.5 


0.272 


0.003 


0.005 


CTEQ6.1 


2155 


123 


109 


578.9 


26.6 


26.8 


0.269 


0.002 


0.007 


MRST2004 (NNLO) 


2193 


41 


45 


598.0 


12.3 


14.5 


0.269 


0.004 


0.003 


MRST2004 (NLO) 


2223 


42 


46 


588.5 


12.1 


14.3 


0.268 


0.004 


0.003 



Table 7: Cross-sections cr, and acceptances A, with asymmetric Hessian uncertainties as calculated 
using several recent PDF sets for the three cut regions defined in Table 0. 



MC@NLO at scale fi^ = = Mz and interfaced to PHOTOS to add final state QED 
radiation. The missing electroweak contribution may then be inferred from HORACE as in 
Sec. ^ For these errors we take those resulting from the tight cut set in Table |l|; these cuts 
are considered the most representative of likely analysis cuts for the LHC experiments. 

QCD uncertainties may be divided into two main classes. If the NNLO K-iactor is set 
to 1, there is a missing NNLO contribution (5nnlo = K — 1. Since K has residual NNLO 
scale dependence, we must also take this into account and write K = 1 + Jnnlo i Jscaie- 
The factor 5nnlo can be inferred from Table |^ and 5 scale can be inferred from half the 
scale variation of the NNLO entries in Table For example, for the total cross-section, we 
can infer that (^nnlo = -(0.98 ± 0.22)%, while (5 scale = (0.64 ± 0.16)%. Both errors in the 
MC@NLO result are of order 1%. The original NLO scale dependence estimate of 5.13% is 
no longer relevant because K is now known, even though it may be set to 1 in a particular 
calculation. Similar estimates can be made for the two sources of error in each of the cuts. 

Both classes of QCD errors are also associated with a "technical precision" due to 
limitations of the computing tools used to evaluate them. Significant improvements in the 
NNLO precision could be obtained if a program with faster convergence were available. We 
therefore include an "error on the error" for the QCD errors, and propagate these through 
in the usual fashion to derive a final accuracy for the total QCD uncertainty estimate. 
This sets a limitation on how much the NLO calculation can realistically be improved 
using currently available NNLO results. 
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QCD Uncertainties (%) 



Cross-Section Ao" 


Uncertainty 


^tot 


Cut 1 


Cut 2 


Cut 3 


Missing NNLO 
Scale Dependence 


-0.98 ±0.22 
0.64 ±0.16 


1.42 ±0.77 
0.51 ±0.54 


-0.98 ±2.80 
1.16 ± 1.76 


-1.33 ±0.62 
0.83 ± 0.46 


Total 


1.17 ±0.20 


1.51 ±0.75 


1.52 ±2.25 


1.57 ±0.58 


Error in Acceptance (AA) 


Uncertainty 




Cut 1 


Cut 2 


Cut 3 


Missing NNLO 
Scale Dependence 




2.43 ±0.81 
0.77 ±0.56 


0.00 ± 2.68 
0.83 ± 1.77 


-0.05 ±0.66 
1.15 ±0.48 


Total 




2.55 ±0.79 


0.83 ± 1.77 


2.55 ±0.79 



Table 8: Summary of QCD uncertainties Act in the cross-sections and A A in the acceptances 
relative to the un-cut cross-section ct*°* for Mu > 40 GeV/c^. The three cuts are described in 
Table ||. Missing NNLO is shown with a sign, because it has been calculated. 

These contributions to QCD errors are summarized in Table |8| for the total cross-section 
and the three cuts of Table |2|. Results are shown both for the three cut cross-sections and 
their ratio for the total cross-section with > 40GeV/c^, and the errors are assumed to 
be uncorrelated. 

If the ii'- factor had not been calculated at NNLO, the error of the NLO cross-sections 
could have been roughly estimated from half the width of the scale-dependence band, or 
half the NLO results for Aa/W in Table |, giving uncertainties of 2.2 — 2.6%. The errors 
calculated from the K factors are in within the limits these expectations, up to the technical 
precision of the calculation. A similar error NLO estimate for the error in the acceptance 
based on Table ^ would predict at most 0.33% missing NNLO. The actual missing NNLO 
in Cut 1 was found to be larger than this, but the other two cuts have very small missing 
NNLO, to within the technical precision of the calculation. 

The final contribution to the total error considered here is the uncertainty from the 
PDFs. This may be extracted from the results of Sec. ^ by taking the errors from the 
CTEQ6.5 results for Cut 1 (see Table ||). The errors are asymmetric, so we take the largest 
of the two (up or down) uncertainties as the total fractional error for the PDF calculation. 
We choose the first cut set, since it is the most representative of likely analysis cuts at the 
LHC experiments. CTEQ errors, rather than the MRST errors, are used because they give 
a more conservative estimate. The difference between the results obtained by the latest 
CTEQ and MRST PDFs is less than the maximum error quoted for CTEQ for all three 
cut regions. 

The errors are added in quadrature, assuming no correlations, and the results are given 
in Table |9|. It should be remarked that there is no concensus on the best way to combine 
the errors in Table |9|, so this total must be considered an estimate. The QCD error is 
taken for Cut 1 for the same reasons as given above. In addition, as we have discussed 
above, we propagate the "error on the error" for each of the contributions in order to 
have some reasonable estimate of the accuracy of the quoted total theoretical uncertainty. 
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The exception to this is the PDF error, which can be considered as an upper hmit on 
the uncertainty. We conclude that the event generator MC@NLO interfaced to PHOTOS 
should be sufficient to guarantee an overall theoretical uncertainty on the Z production 
due to higher order calculation, PDFs, and renormalization scale at the level of 4% for the 
total cross-section and at approximately 3% for the acceptance. 

Total Theoretical Uncertainty (%) 



Uncertainty 


Cross-Section Au 


Acceptance Ayl 


Missing 0(a) EWK 
Total QCD Uncertainty 
PDF Uncertainty 


0.38 ±0.26 
1.51 ±0.75 
3.79 


0.96 ±0.21 
2.55 ±0.79 
1.32 


Total Uncertainty 


4.1 ±0.3 


3.0 ±0.7 



Table 9: Total theoretical uncertainty on the Z production cross-section Act, and acceptances A^. 



Z production will provide a valuable tool for studying QCD, measuring precision elec- 
troweak physics, and monitoring the luminosity. As the luminosity increases, the large 
statistics will permit a further improvement in the systematic uncertainties due to the 
PDFs. Adding complete 0(a) EWK corrections to the event generator would eliminate 
most of the EWK uncertainty, and incorporating NNLO QCD corrections would substan- 
tially reduce the QCD uncertainties. 

Reaching a combined precision of 1%, as desired in the later stages of analysis at high 
integrated luminosity, will require new tools. In addition to improved PDFs, an event gen- 
erator combining NNLO QCD with complete 0(a) EWK corrections will be needed, with 
exponentiation in appropriate regimes, and adequate convergence properties to technically 



reach the required precision. 1 27] Measurement strategies have been proposed that may mit- 
igate some of the effects of systematic errors on the precision of the Z measurements. A 
combination of improved calculations and improved meaurements will be needed to permit 
the desired precision to be reached as the integrated luminosity increases to a point where 
it is needed. 
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(c) 



Figure 11: Comparison of Z/j* £^£ {£ — e or /i) acceptances A, with several recent PDF 
calculations for acceptance regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table |. The 
left-hand plots show the total acceptance and the right hand plots show the fractional error on the 
acceptance. 
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Figure 12: The Z/'^* £^£~ cross-section cr = e or /i), as a function of the p^, cut for 
acceptance regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table |[ For each acceptance 
region we fix the invariant mass and \r]\ cuts at their specified values, and vary only the cut. 
The figures on the right show the relative errors in the cross sections. 
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Figure 13: The Zj^* t^t~ cross-section cr (£ = e or /i), as a function of the \'r]\ cut for 
acceptance regions (a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table |[ For each acceptance 
region we fix the invariant mass and p^cuts at their specified values, and vary only the \r]\ cut. The 
figures on the right show the relative errors in the cross sections. 
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Figure 14: The Zj^* t^t~ acceptances A, as a function of the -p^ cut for acceptance regions 
(a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table ^. For each acceptance region we fix the 
invariant mass and |?7| cuts at their specified values, and vary only the cut. The figures on the 
right show the relative errors in the acceptances. 
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(c) 



Figure 15: The Z/'y* t^(,~ acceptances A, as a function of tiie \r]\ cut for acceptance regions 
(a) Cut 1, (b) Cut 2, and (c) Cut 3, as defined in Table ^. For eacli acceptance region we fix tlie 
invariant mass and p^^cuts at their specified values, and vary only the |ry| cut. The figures on the 
right show the relative errors in the acceptances. 
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